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Abstract 

The i?-matrix formalism for the construction of integrable systems with infinitely 
many degrees of freedom is reviewed. Its application to Poisson, noncommutative and 
loop algebras as well as central extension procedure are presented. The theory is de- 
veloped for (1 -|- l)-dimensional case where the space variable belongs either to M or 
to various discrete sets. Then, the extension onto (2 -|- l)-dimensional case is made, 
when the second space variable belongs to M. The formalism presented contains many 
proofs and important details to make it self-contained and complete. The general the- 
ory is applied to several infinite dimensional Lie algebras in order to construct both 
dispersionless and dispersive (soliton) integrable field systems. 
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1 Introduction 

One of the most characteristic features of integrable nonlinear systems is the existence of 
bi-Hamiltonian structures. This ingenious concept was introduced by F. Magri [52] in 1978. 
From the geometrical point of view, this means that there exists a pair of compatible Poisson 
tensors which allow, using a recursion chain, to generate infinite (in the infinite-dimensional 
case) hierarchies of commuting symmetries and constants of motion being in involution with 
respect to the above Poisson tensors. In order to stress the importance of the bi-Hamiltonian 
structures for evolution systems let us quote L.A. DickejS 

" The existence of two compatible Poisson (or Hamiltonian) structures is a remark- 
able feature of the most, if not all, integrable systems, sometimes it is considered 
as the essence of the integrability." 

Finding a systematic method for construction of integrable systems is one of the most 
important issues in the theory of evolutionary systems. It is well known that a very powerful 
tool called the classical R-matrix formalism [Tlj proved to be very fruitful in the systematic 
construction of the field and lattice soliton systems as well as dispersionless systems. The 
crucial point of the formalism is the observation that integrable dynamical systems can be 
obtained from the Lax equations on appropriate Lie algebras. The greatest advantage of 
this formalism, besides the possibility of systematic construction of the integrable systems, 
is the construction of bi-Hamiltonian structures and (infinite) hierarchies of symmetries and 
conserved quantities. 

The goal of this article is to present the formalism of classical i?-matrices applied to infinite- 
dimensional Lie algebras in order to construct integrable systems with infinitely many degrees 
of freedom and related tensor invariants like conserved functions, symmetries, Poisson tensors, 
recursion operators and so on. 

In the first part of the paper we present in a systematic fashion the concept of i?-matrix 
formalism with many proofs and important details to make the text self-contained and com- 
plete. First of all we present the basics of the formalism and then illustrate its approach to 
the Poisson and noncommutative algebras, with the aim of construction of integrable hier- 
archies with multi-Hamiltonian structures. This theory is developed for (l-l-l)-dimensional 
systems, where one variable is the evolution parameter t (time), which in presented construc- 
tion always belongs to M, and the spatial variable, which may belong to K = M or C but also 
to various sets of discrete numbers, so we cover the cases like continuous dynamics, lattice 
dynamics and g-dynamics in a single formalism. Then, we show how to extend the whole 
formalism over additional spatial variable via the so-called central extension approach. The 
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advantage of this construction is the possibihty, in principle, of finding a (2+l)-dimensional 
counterpart for an arbitrary (l+l)-dimensional system derived via the i?-matrix formahsm, 
with the only limitation that the additional space coordinate is always continuous (belongs 
to M or S^) rather than discrete. Finally we apply the i?-matrix formalism to loop algebras 
together with central extension procedure. 

In the second part of the paper we apply the formalism developed in the first part to 
several important Lie algebras. Some of them allow to construct dispersionless integrable 
systems while the others lead to dispersive (soliton) integrable systems. In the second case 
we construct integrable systems with continuous number of degrees of freedom as well as 
those with discrete number of degrees of freedom. As each algebra leads to a huge number of 
specific integrable systems, in this paper we present explicitly only few of them, referring the 
reader to the cited references for further examples. 

In the following review classical i?-matrix theory is presented in the framework of an 
appropriate Lie algebras, with the aim of construction of integrable bi-hamiltonian systems. 
However, the formalism is significantly more powerful as it is intimately connected with 
factorization and Riemann-Hilbert problem for the related Lie groups 

2 Classical i?-matrix theory 

In this chapter we will present a unified approach to the construction of integrable evolu- 
tion equations together with their (multi-)Hamiltonian structures. The idea originates from 
the pioneering article [3lj by Gelfand and Dickey, where they presented a construction of 
Hamiltonian soliton systems of KdV type using pseudo-differential operators. Next, Adler [2] 
showed how to construct the bi-Hamiltonian structures for the above soliton systems using 
the method based on the Kostant-Symes theorem obtained independently in [H] and |83j . 
Later the abstract formalism of classical /^-matrices appeared in the paper [77] by Sklyanin as 
an intermediate step within the inverse quantum scattering method. The present version of 
the formalism was given in [711 [71] . In [SOI [22] it was shown that there are in fact three natu- 
ral Poisson brackets associated with classical i?-structures. Quite recently Li [lU] considered 
the classical i?-matrix theory on the so-called (commutative) Poisson algebras. This approach 
leads to the construction of multi-Hamiltonian systems of hydro dynamic (dispersionless) type. 

2.1 Classical i?- matrices 

Let g be a Lie algebra over the field IK of complex or real numbers, IK = C or M, that is, g 
is equipped with a bilinear operation [-, ■] : x g g, called a Lie bracket, which is skew- 
symmetric and satisfies the Jacobi identity. The Lie bracket [■, ■] defines the adjoint action of 
g on g: adab = [a,b]. 

Definition 2.1 A linear map i? : g — g such that the operation 

[a, b]j^ := [Ra, b] + [a, Kb] a,b e Q (2.1) 
defines another Lie bracket on g is called the classical R-matrix. 
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The skew-symmetry of (12.1 1) is obvious. As for the Jacobi identity for (12. ip . we find that 

= [a, [b, c]^]^ + c.p. = [Ra, [Rb, c]] + [Ra, [b, Rc]] + [a, R [b, c] J + c.p. 
= [Rb, [Rc, a]] + [Rc, [a, Rb]] + [a, R [b, c] J + c.p. 

= [a,R[b,c]j^-[Rb,Rc]] + c.p., (2.2) 

where c.p. stands for cychc permutations within the triple {a, b, c} G g, and the last equality 
follows from the Jacobi identity for [■,■]. Hence, a sufficient condition for i? to be a classical 
i?- matrix is to satisfy the so-called (modified) Yang-Baxter equation, YB{a): 

[Ra,Rb]- R[a,b]j^ + a[a,b] = 0, (2.3) 

where a is a number from K. There are only two relevant cases of YB(a), namely a ^ and 
a = 0, as all Yang-Baxter equations with a ^ are equivalent up to a rescaling of a. 

Definition 2.2 A linear operator A : g ^ g is called intertwining if Ao ada = ad^ o A, i.e., 
if A [a, b] = [Aa, b] = [a, Ab] for any a,b E g. 

Proposition 2.3 J7^ If R is a classical R-matrix and A an intertwining operator, then RoA 
also is a classical R-matrix. 



Proof. We have 



RA [a, 6]^^ = RA [RAa, b] + RA [a, RAb] 

= R [RAa, Ab] + R [Aa, RAb] = R [Aa, Ab]j^ . 

Hence, 

[RAa, RAb] - RA [a, 6]^^ = [RAa, RAb] - R [Aa, Ab]^^ 

and the Jacobi identity for [-, ■]ra with respect to the elements a,b,c E g reduces to the Jacobi 
identity for [■,-]r with respect to the elements Aa, Ab, Ac, see (12.21) . □ 

Notice that a linear combination of intertwining operators again is an intertwining oper- 
ator. 



2.2 Lax hierarchy 

In this section we present the classical i?-matrix formalism for the class of Lie algebras for 
which the Lie bracket additionally satisfies the Leibniz rule. Later, while considering the loop 
algebras in Section 12.101 we shall drop this extra condition. 

Assume that the Lie algebra g is also an algebra with respect to an associative multipli- 
cation such that 

ada(6c) = ada(6)c + 6ada(c) <^==^ [a,bc] = [a,b]c + b[a,c] (2.4) 

the Leibniz rule holds, i.e., the Lie bracket [■, ■] is a derivation with respect to the multiplica- 
tion. Notice that this condition is satisfied automatically in the case of a commutative algebra 
g when the Lie bracket is given by a finite-dimensional Poisson bracket, as well as in the case 
of a non- commutative algebra g with the Lie bracket given by the commutator. 
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As a consequence of the above assumption, any well-defined smooth mapping X : q ^ q, 
L I— > X{L), is an invariant of the Lie bracket, i.e., 

[X{L),L]=0. (2.5) 

The smoothness of the mapping X means that its differential and directional derivatives exist 
and are well-defined. 

Proposition 2.4 Smooth invariant functions Xn{L) generate a hierarchy of vector fields on 
g of the form 

Lt^ = [RXn{L),L] Leg n e Z, (2.6) 

where t„, are evolution parameters. Assuming that a classical R-matrix R commutes with the 
directional derivatives with respect to all ( \2.6\} . the Yang- Baxter equation (12.31) is a sufficient 
condition for the pairwise commutativity of the vector fields (12.61) . 

Proof. The directional derivative of a smooth function F : g ^ g in the direction of (12.61) 
is given by 

F{L),„ = F{L)' [L,J = [RX^{L), F{L)] , (2.7) 
which follows from the Leibniz rule (12. 4p . Thus, one finds that 

(L,J,^ - {L,J,^ = [RXUL),L],^ - [RX^{L),L],^ 

= [{RXUL)),^ - {RX^{L)\^,L] + [RX^{L), [i?X„,(L), L]] - [i?X„(L), [RX^{L),L\] 
= [{RX^{L)),^ - (RX^iL)),^ + [RXUL), RX^{L)] , L] . (2.8) 

Hence, for the pairwise commutativity of vector fields (12.61) it suffices that the so-called zero- 
curvature equations 

{RXr^iL)),^ - (RX^iL)),^ + [RXUL), RX^{L)] = (2.9) 

hold. 

The assumption that R commutes with directional derivatives implies that it commutes 
with the derivatives with respect to evolution parameters, i.e., 

{RL\^ = RU^. (2.10) 

Thus, the right hand side of (12. 9p becomes 

R{XM\. -R{X^{L)\^^ + [RXUL),RX^{L)] = 

'^1^ i?[i?X„(L),X„(L)] -R[RXUL),X4L)] + [RX^L), RX^L)] 

= [RXrn{L),RX4L)]-R[Xm.iL),X4L)]^. (2.11) 

Now, if the i?-matrix R satisfies the Yang-Baxter equation (12.31) then the last expression is 
equal to —a [Xm{L), Xn{L)] = 0, and the result follows. □ 

The hierarchy (12. 6p is called the Lax hierarchy and L is called the Lax operator or the Lax 
function depending on the nature of a given Lie algebra g. Notice that the assumption that R 
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commutes with directional derivatives is an important condition although is not enunciated 
explicitly in most works on the i?-matrices. 

The natural choice for invariant smooth functions are the power functions X„(L) = 
that are always well-defined on g. One can consider less trivial functions, for example the 
logarithmic ones, like X{L) = InL, but only when they have proper interpretation in g. 

It is natural to ask when the abstract Lax hierarchy (12.61) represents a "real" hierarchy of 
integrable evolution systems on a suitable function space constituting an infinite-dimensional 
phase space L(. This occurs when we can construct an embedding map l : U ^ q, which 
induces the differential structure on q. Notice that for l being an embedding its differential 
t' : V — is an injective map, where V is a space of vector fields on U. In such a case the Lax 
hierarchy fl2.6l) can be pulled back to the original function space by The symmetries from 
the Lax hierarchy fl2.6p represent compatible evolution systems when the left- and right-hand 
sides of (12.61) span the same subspace of g. So, the Lax element L of g has to be chosen in a 
suitable fashion. 



2.3 Simplest i?-matrices 

The simplest way to obtain a classical i?-matrix is to decompose a given Lie algebra into Lie 
subalgebras. Thus, assume that a Lie algebra g can be split into a (vector) direct sum of Lie 
subalgebras g+ and g_, i.e., 

S = 0+©S- [0±,0±]cg± g+ng_ = 0. 

It is important to stress that we do not require that [g+,g_] = 0. 

Upon denoting the projections onto the subalgebras in question by P±, we define a linear 
map i? : g ^ g as 

i?=l(P+-P_). (2.12) 
Using the equality P+ + P_ = 1 (12.121) can be represented in the following equivalent forms: 

P = P+-i = i-P_. (2.13) 

Let a± := P±(a) for a G g. Then 

1 1 

[a,6]^ = [a+,6+] - [a_,6_] =^ P [a, 6]^ = - [a+, 6+] + - [a_, 6_] 

and 1111 

[Ra, Rb] = - [a+, b+] - - [a+, 6_] - ^[a~, b+] + ^ b^] . 

Hence, the map (12.121) satisfies the Yang-Baxter equation (12.31) for a = ^ and is a well-defined 
classical P-matrix. This is the simplest and the most common example of a well-defined 
P-matrix. 

For instance, the Lax hierarchy (12. 6p for the P-matrix (I2.12p . following from the decom- 
position of a Lie algebra into Lie subalgebras, takes the form 

L,„ = [(X„(L))^,L]=-[(X„(L))_,L]. 

It can be written in two equivalent ways because (12.130 holds. 
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The construction of majority of known integrable systems within the formahsm presented 
above is based on the classical i?-matrices that follow from the double decomposition 02.121) 
of Lie algebras into Lie subalgebras. In [571 ESI EHl E5] the authors considered deformations 
of (12.121) that preserve the Yang-Baxter equation and originate from a triple decomposition 
of a given Lie algebra; see also [87] for multiple decompositions of Lie algebras. 



2.4 Lie-Poisson structures 

Let Q* be a (regular) dual of a given Lie algebra q and (■,■): g* x g —> IK be the usual duality 
pairing. The co-adjoint action ad* of g on g* is defined through the relation 

(ad:r/, b) + (r/, ad,6) = ^ (ad>, b) = - (r/, [a, b]) , (2.14) 

where a, & G g and 77 G g* 

Let this time l : U ^ g* he the embedding of the original phase space into the dual Lie 
algebra. Then every functional F : U can be extended to a smooth function on g*. 

Therefore, let C°°(g*) be the space of all smooth functions on g* of the form F o : g* — s- K, 
where F e C°°{U). Then the differentials dF{'q) of F{ri) e C°°{q*) at the point r] e g* belong 
to g as they can be evaluated using the relation 



de 



= {^,dFiv)) eeg* eeK. (2.15) 

e=0 



Moreover, the form of differentials dF G g* has to be such that the duality pairing between 
g and its dual g* coincides with the duality map between vector fields and one-forms on the 
original infinite-dimensional function phase space U. Indeed, 



iVudF) 



Y,T-iui)tdx, (2.16) 

i=0 * 



where rjt G g* is a vector field on g*, F{r]) G C°°{g*) is a functional depending on the dynamical 
fields Ui from the phase space U, is the variational derivative of F with respect to field 
variable Ui and in the case of discrete functionals the symbol of integration in f l2.16p must be 
replaced by an appropriate summation. We also have the relation 

(e,rfF'M) = (r/,rfF'[e]), (2.17) 

which is equivalent to the vanishing of the square of the exterior differential, i.e., d'^F = 0. 

We also make an additional assumption that the Lie bracket in g is such that directional 
derivative along an arbitrary ^ G g* is a derivation of the Lie bracket. This means that the 
following relation holds: 

[a,bn^] = [aUb] + [a,b'[^]]. (2.18) 

Theorem 2.5 There exists a Poisson bracket on the space of smooth functions on a dual 
algebra g* , which is induced by the Lie bracket on g. This Poisson bracket is defined as 
follows: 

{H,F}{ri):={rj,[dF,dH]) r, e g* H,FeC^{g*). (2.19) 
Lemma 2.6 The differential of (12.191) is given by 

d {H, F} = [dF, dH] + dF' [ad^^r/] - dH' [ad^^r/] . (2.20) 
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Proof. By f lZT^ we find that 

{H, F}' [e] = (ry'le], [dF, dH]) + (r/, rf//] + [dF, dH'[i]]) 

'^i® (e, + {^d*,Hr,,dF%]) - {^d*,,r,,dH'[i]) 

and the result follows. □ 



Proof of Theorem 12.51 Bilinearity and skew-symmetry of (12.191) are obvious, so we only 
have to prove the Jacobi identity: 

{F, {G, if}} + c.p. = {r], [d {G, H} , dH] + c.p.) 



(r/, [[rfii, dG] ,dF]) + (ad^^r/, (iii' [ad^^r/]) - (ad^^r/, dG' [ad^^r/]) + c.p. 

I 



•"^iz-P- (77, [[dH, dG] ,dF]) + (ad;j^?7, rfii' [ad^c^]) - {^(KigV, dH' [ad;^^.??]) + c.p. 



[t], [[dH,dG],dF]+c.p.) = 0, 
where the last equality follows from the Jacobi identity for [-,■]. □ 

The bracket (12.191) is called a (natural) Lie-Poisson bracket and was originally discovered 
by Sophus Lie. Its modern formulation is due to Berezin [6] as well as Kirillov and Kostant |41j . 

Now assume that we have an additional Lie bracket (12. ip on g defined through the classical 
i?-matrix such that (I2.10p is valid. Then (12. ip satisfies the condition (I2.18p . As a result, there 
is another well-defined (by Theorem 12. 5p Lie-Poisson bracket on the space of scalar fields 
C°°(0*): 

{H,F}^{r^):={r^,[dF,dH]^) rj E g* H^FeC^iQ*). (2.21) 

Using (I2.14P we find that the associated Poisson operators at 77 G g*, the one for the 
natural Lie-Poisson bracket (I2.19P and the second one, (I2.2ip . have the form 

{H, F} = {ndH, dF) ^ tt : dH ^ ad*aHV 

{H, F}^ = {iTRdH, dF) ^ TlR-.dH^ ^^RdHV + R*^^*dHV. 

where the adjoint of R is defined by the relation (-R*//, a) = (77, Ra), where 77 G g* and a G g. 
The following theorem constitutes the essence of the classical i?-matrix formalism. 

Theorem 2.7 (| I74] ) The Casimir functions Cn G C°°{q*) of the natural Lie-Poisson bracket 
(I2.19P are in involution with respect to the Lie-Poisson bracket (I2.2ip induced by (12.10 . More- 
over, Cn generate a hierarchy of vector fields on g* : 

Vt„ = T^RdCniv) = ad*j^dc„v V e g*. (2.22) 

The evolution systems (12.220 pairwise commute, i.e., {rit^)t„ = {Vt„)tm' ^'^^ ^'^^ Hamiltonian 
with respect to (I2.2ip . Moreover, any equation from (I2.22p admits all Casimir functions Cn 
of Ii2.19\) as integrals of motion. 
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Proof. The Casimir functions C„ of the natural Lie-Poisson bracket fl2.19p satisfy the fol- 
lowing condition 

VFgC~(0*) {F,Cn} = 

that is, their differentials are ad*-invariant. Hence, they are in involution with respect to the 
Lie-Poisson bracket f l2.2ip . i.e., {C„, Cm}^ = 0. Now, as Hfid is a Lie algebra homomorphism 
from the Poisson algebra of smooth functions on q* to the Lie algebra of vector fields on g*, 
commutativity of Hamiltonian vector fields fl2.22p . with the Casimir functions as Hamiltoni- 
ans, readily follows. □ 

In fact, when the i?-matrix follows from the decomposition of an Lie algebra into sum of 
Lie subalgebras, i.e. is given by (12.121) . Theorem 12.71 can be considered as a generalization of 
Kostant-Symes theorem [HI |H3] • 

The construction of Casimir functions C„ and related dynamical systems (I2.22p on the 
dual Lie algebra g* is, in contrast with (12. 6p . often inconvenient and impractical. Thus, a 
formulation of a similar theory on g instead on g* is often justified. This can be done when 
one can identify g* with g by means of a suitable scalar product. 



2.5 74(i-invariant scalar products 

We restrict our further considerations to the Lie algebras g for which their duals g* can be 
identified with g through a duality map. So, we assume the existence of a bi-linear scalar 
product 

(■,-)g:Sxg^K (2.23) 

on g, and we assume this product to be symmetric: (a, b)^= {b, a)^, and non-degenerate, that 
is, a = is the only element of g that satisfies (a, b)^ = for all & G g. Then we can identify 
g* with g (g* = g) by setting (77, b) = (c, b)^, V6 G g, where 77 G g* is identified with c G g. We 
also make an additional assumption that the symmetric product (I2.23P is ad-invariant, i.e., 

{[a,c],b)^ + {c,[a,b])^ = 0. (2.24) 

This is a counterpart of the relation (I2.14p . Thus, if 77 G g* is identified with c G g we have 
(ad*r;, b) = ([a, c] ,b)^ and one identifies ad*// G g* with ad^c G g. 

In fact, by virtue of the scheme presented in the previous section all equations from the 
hierarchy (12. 6p are Hamiltonian. Since g* = g, the Lie-Poisson brackets (I2.19P and (I2.2ip on 
the space of scalar fields C°°(g = g*) at L G g take the form 

{H, F} = (L, [dF, dH])^ = {dF, TidH)^ ^ ndH = [dH, L] 

{H, F}^ = (L, [dF, dH]^)^ = {dF, irndH)^ ^ undH = [RdH, L] + R* [dH, L] , (2.25) 

where now R* is defined by the relation {R*a, &)g = (a, Rb)^. 

Differentials of the Casimir functions Cn{L) G C°°(g) of the natural Lie-Poisson bracket 
are invariants of the Lie bracket, i.e., [dCn{L), L] = 0. Obviously, the Casimir functions are 
still in involution with respect to the second Lie-Poisson bracket defined by R and generate 
pairwise commuting Hamiltonian vector fields of the form 



Lt^=7iRdC^{L) = [RdCn,L] 



(2.26) 
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Notice that the Lax hierarchy fl2.6l) coincides with fl2.26p for Xn = dCn- Moreover, it 
follows that if there exists a symmetric, non-degenerate and ad-invariant product on g then 
the Yang-Baxter equation f l2.3p is not a necessary condition for the commutativity of vector 
fields from the Lax hierarchy (12. 6p . However, if (12.30 is not satisfied then the zero-curvature 
equations (12. 9p will not be automatically satisfied as well. 

The simplest way to define an appropriate scalar product on a Lie algebra g is to use a 
trace form Tr : g — > K such that the scalar product 

(a,6)g := Tr(a6) a, 6 G (2.27) 

is nondegenerate. In this case the symmetry of (12.270 entails that 

Tr {ah) = Tr (ba) . (2.28) 



Lemma 2.8 Let Tr : g ^ IK &e a trace form defining a symmetric and nondegenerate scalar 
product (12.271) such that the trace of Lie bracket vanishes, i.e. Tr [a, 6] = for all a,b E g. 
Then the condition (12.41) for the Lie bracket to be a derivation with respect to the multiplication 
is a sufficient condition for (12.270 to be a.d-invariant. 

Moreover, if the Lie bracket in g is given by the commutator, [a, b] = ab — ba, then the 
ad-invariance follows from the associativity of the multiplication in g. 

Proof. The first part of the lemma is immediate, as 

([a, c] , 6)g + (c, [a, b])^ = Ti {[a, c]b + c [a, b]) = Tr [a, cb] = 0, 

where we used the assumptions from the proposition. The second statement of the lemma 
follows immediately from the assumption and (12.280 . □ 

Under the assumption that we have an algebra g such that a Lie bracket is a derivation of a 
multiplication (12. 4p and g is endowed with a trace form inducing a nondegenerate ad-invariant 
scalar product (I2.27p . the most natural Casimir functions C„(L) G C°°(g) of the Lie-Poisson 
bracket (I2.19P are given by the traces of powers of L, i.e., 

Cn{L) = ^^Tr ^ dCn = L'^ n ^ -I. (2.29) 

The associated differentials are found from the expression (12.150 . which can be now reduced 
to 

^F(L) = (L„ dF)^ = Tr (L, dF) L G g, (2.30) 
where t is an evolution parameter associated with a vector field Lt on g. 

2.6 Hamiltonian structures on Poisson algebras 

Definition 2.9 Let A be a commutative, associative algebra with unit 1. // there is a Lie 
bracket on A such that for each element a E A the operator ad^ : b i— > {a, b} is a derivation 
of the multiplication, i.e. {a, be] = {a, b}c + b{a, c} then {A, {■,■}) is called a Poisson algebra 
and the bracket {■, ■} is a Poisson bracket. 
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Thus, the Poisson algebras are Lie algebras with the Lie bracket [■, ■] := {■, ■} endowed with 
an additional structure. Of course we should not confuse the above bracket with the Poisson 
brackets in the algebra of scalar fields. It will follow easily from the context which bracket is 
used. In the case of the Poisson algebra A a classical i?-matrix defines the second Lie product 
on A but not the Poisson bracket; in general this would not be possible. 

Theorem 2.10 Let A be a Poisson algebra with the Poisson bracket {■,■} and a non- 
degenerate a.d-invariant scalar product (■, ■)^ such that the operation of multiplication is sym- 
metric with respect to the latter, i.e., {ah, c)^ = (a, 6c)^ for all a,b,c & A. Assume that R is 
a classical R-matrix such that f l2.10p holds. 
Then for any integer n ^ the formula 

{H, F}„ = (L, {RiVdF), dH} + {dF, R{U'dH)})j^ , (2.31) 

where H, F are smooth functions on A, defines a Poisson structure on A. Moreover, all 
brackets {■,-}n o,re compatible. 

An important property that classical i?-matrices commute with differentials of smooth maps 
from ^ to or equivalently satisfy (12.101) , is used in the proof of Theorem 4.2 of [19] , although 
it is not explicitly stated there. In fact, the existence of scalar product being symmetric with 
respect to the multiplication, {ab, c) = (a, be), entails existence of a trace form on A. Setting 
c = 1 we have {ab, 1) = (a, b). Thus, the trace can be defined as Tr(a) := (a, 1) = (1, a). 

The Poisson operators vr^ related to the Poisson brackets (12.311) such that {H, F}^ = 
{dF, TindH) , are given by the following Poisson maps 

TT^: dH ^ {R{U'dH),L} + L^'R* {{dH,L}) n ^ 0. (2.32) 

Notice that the bracket (12.311) with = is just the Lie- Poisson bracket with respect to the 
second Lie bracket on A defined by a classical i?-matrix. Referring to the dependence on L, 
the Poisson maps (I2.32p are called linear for n = 0, quadratic for n = 1, and cubic for n = 2, 
respectively. The Casimir functions C{L) of the natural Lie- Poisson bracket are in involution 
with respect to all Poisson brackets (I2.32p and generate pairwise commuting Hamiltonian 
vector fields of the form 

Lt = -KndC = {R (L'^dC) , L} L e A. 

Taking the most natural Casimir functions (12.290 . defined by traces of powers of L, for 
the Hamiltonians, one finds a hierarchy of evolution equations which are multi-Hamiltonian 
dynamical systems: 

Lt^ = {RdCn, L} = TiodCn = TiidCn-i = ■■■ = nidCn-i = . • • , (2.33) 

where C„ are such that dCn = LdCn-i- For any /?- matrix any two evolution equations in the 
hierarchy (I2.33P commute because of involutivity of the Casimir functions Cg. Each equation 
admits all the Casimir functions as conserved quantities in involution. In this sense we will 
consider (I2.33P as a hierarchy of integrable evolution equations. The most natural choice for 
the Casimir functions are the traces of the power functions (12.290 . 
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2.7 Hamiltonian structures on noncommutative algebras 

In this section, in contrast with the previous one, we will consider a noncommutative as- 
sociative algebra g, with unity, for which the Lie structure is defined as a commutator, i.e., 
[a, b] := ab — ba, where a,b E Q. Such a Lie bracket automatically satisfies the required Leibniz 
rule 02.41) . We further assume existence of nondegenerate, symmetric and ad-invariant scalar 
product on g. Let i? be a classical i?-matrix such that fl2.10l) is satisfied. 

In this case the situation is more involved and only three explicit forms of Poisson brackets 
on the space of smooth functions defined by related Poisson tensors are known from 

the literature: 

{H,F}^ = {dF,7r^dH)^ n = 0,1,2. 

These Poisson brackets (or associated tensors) are called linear, quadratic and cubic brackets 
(resp. tensors) for n = 0, 1, 2, respectively. 

The linear one is simply the Lie- Poisson bracket, with respect to the second Lie structure 
on g defined by classical i?-matrix, with the Poisson tensor (12.251) 

TTodH = [RdH, L] + R* [dH, L] , (2.34) 

for which there is no need for additional assumptions. 

In our further considerations we have to assume that the scalar product is symmetric with 
respect to the operation of multiplication, {ab,c) = [a, be). Note that this property implies 
that the scalar product is automatically ad-invariant with respect to the Lie bracket defined 
by the commutator (Lemma 12. 8p . 

The quadratic case is more delicate. A quadratic tensor [82] 

TTidH = Ai{LdH)L - LAiidHL) + S{dHL)L - LS*{LdH) (2.35) 

defines a Poisson tensor if the linear maps ^41^2 : ^ are skew-symmetric, = —^1,2, 
satisfy YB(a) (12.31) for a 7^ and the linear map S : q ^ g with its adjoint S* satisfy 

S ([^a, b] + [a, A2b]) = [Sa, Sb] , S* {[A^a, b] + [a, A^b]) = [S*a, S*b] . (2.36) 

In the special case when 

R:=l{R-R*) (2.37) 

satisfies the YB(a), for the same a as R, under the substitution A\ = A2 = R — R* and 
S = S* = R + R* the quadratic Poisson operator fl2.35p reduces to 



-KidH = [R [dH, L]^ , L] + LR* [dH, L] + R* {[dH, L]) L, (2.38) 

where [a, 6]^ := ab + ba and the conditions (12.361) are equivalent to YB(a) for R and R. In 
particular, when R* = —R, these conditions are automatically satisfied as in this case R = R. 

Another special case occurs when the maps Ai^2 and S originate from the decomposition 
of a given classical i?-matrix satisfying YB(a) for a 7^ 



R=liA, + S) = l{A2 + Sn, 
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where 2 are skew-symmetric. Then, the conditions fl2.36p imply that both Ai and A2 
satisfy YB[a) for the same value of a as R, [STj. Hence, in this case we only have to check 
the conditions f l2.36p for (12.351) to be a Poisson operator. The latter now takes the form 

TTi-.dH^ 2R (LdH) L - 2LR (dHL) + S {[dH, L]) L + LS* [dH, L] . (2.39) 

Finally, the cubic tensor 112 takes the simple form [59] 

TT2:dH^ [R (LdHL) , L] + LR* {[dH, L]) L 

and is Poisson without any further additional assumptions. 

Once again, taking the Casimir functions (I2.29P defined by the traces of powers of L for the 
Hamiltonians yields a hierarchy of evolution equations which are tri-Hamiltonian dynamical 
systems, 

Lt„ = [RdCn, L] = nodCn = -KldCn^i = TT2dCn-2, (2.40) 

where C„ are such that dCn = LdCn-i- We assumed that tt2 in (12.401) is given by (I2.38P or 
(12.390 . In the first case all three Poisson tensors in (12.400 are automatically compatible. In 
the second case this has to be checked separately. 

2.8 Central Extension Approach 

Let be a Lie algebra with the Lie bracket [■,■]. Consider its extension g := g © IK with the 
Lie bracket given by 

[(a,a),(6,/5)] = ad(a,«)(6,/5) := ([a,6],tu(a,6)) a, 6 G g a,/3GK, (2.41) 

where ad is the associated adjoint action. It is readily seen that (12.410 is a well-defined Lie 
bracket if and only if u; is a two- co cycle. 

Definition 2.11 A two-cocycle on g is a bilinear map : g x g — > K such that 
(i) it is skew- symmetric: uj (a, h) = —uj [b, a), 
(a) and it satisfies the cyclic condition: 

uj{[a,b],c) +uj{[c,a],b) +uj{[b,c\,a) = 0, (2.42) 

where a,b,c & g. 

Notice that (0, a) G g commute with respect to (I2.4ip with all other elements from g and 
hence g can be identified with g © a for fixed a G K. The value a is often called a charge. In 
fact g = g/c, where c = {(0, a) G g : a G K} is in the center of g, and thus the Lie algebra g 
is called a central extension of g. 

Assume now (for simplicity) that g can be identified with g* through a non-degenerate 
symmetric scalar product (12.230 . Then this product can be extended to the algebra g in the 
following fashion: 



{{a,a),{b,P)V:={a,b) +aP a, 6 G g a,PeK. 



(2.43) 
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Of course, fl2.43p is symmetric and an important fact is that it preserves non-degeneracy. 
Thus, Q* can be identified with g. 

Usually two-cocycles are defined through the scalar product on g and a linear map : 
g — s> g such that 

u;(a,6) = (a,0(6))g. (2.44) 
If the linear map : g — > g is skew-symmetric, i.e., (p* = — </>, and the following condition 

([a, 6]) = adX6) - adl(j){a) (2.45) 

is satisfied, then it is called a one-cocycle. 

Proposition 2.12 The bilinear form given by f l2.44p is a two-cocycle if and only if (p is a 
one-cocycle. Moreover, if the symmetric product on g is ad-invariant f l2.24p then f l2.44p is a 
two-cocycle if and only if the skew- symmetric (p is a derivation of the Lie bracket [-, ■] in g, 
i.e., 

0( [a, 6]) = [0(a), 6] + [a, 0(6)] (2.46) 

holds. 



Proof. It is clear that fl2.44p is skew-symmetric if and only if 0* = —0. The cyclic condition 
(12.421) for skew- symmetric has the form 

uj {[a, b],c) + c.p. = - (0 {[a, b]) , c)g + {[b, c] , 0(a))g + ([c, a] , 0(6))^ 

= - (0 ([a, b]) , c)g - (ad:0(a), c)^ + (ad:0(6), c)^ = 0, 

where we used the definition (12.141) of coadjoint action. Now, since the symmetric product 
on g is non-degenerate, the cyclic condition is equivalent to (I2.45p . For the ad-invariant sym- 
metric product (I2.46P follows from (12.450 since ad* is identified with ad. □ 

The adjoint action (I2.4ip does not depend on a, thus in fact it defines adjoint action of 
g on g. Hence, we can omit dependence on the charge and write ad^ = ad^a,a)- When a 
given two-cocycle is defined by means of a one-cocycle, i.e., (I2.44p holds, then we can write 
explicitly the coadjoint action ad , since 

(adf,{a, a), (c, 7) := - (^(a, a), adfe(c, 7) = - (a, adbc)^ - a {b, 0(c))g 
= (adja + a0(6),c)g = ((ad^a a0(6), 0), (c, 7))^ . 

Hence, we can restrict ad to g and define the central extension of coadjoint action of g* on 
g, i.e., ad^ := ad^ -|- a0(fe) for 6 G g, where the charge a G K is now treated as a parameter. 

According to Theorem 12.51 the Lie bracket (I2.4ip on g defines a Lie-Poisson bracket on the 
space of smooth functions on g* = g, i.e., on C°°(g). This Poisson bracket can be restricted 
to C°°(g) considered as a subspace of C°°(g). Hence, at a point (L, a) G g we have 

{H, F} (L) := ((L, a), [{dF, 0), {dH, 0)])^ = (L, [dF, dH])^ + au {dF, dH) , (2.47) 

where H,F & C°°{9) and a G K. The Poisson bracket (I2.47P is a central extension of the 
natural Lie-Poisson bracket on C°°(g) generated by the Lie algebra structure on g. When 
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the two-cocycle is given by fl2.44p . then the associated Poisson tensor ir such that {H,F} = 
{dF, irdH)^ has the form iidH = ad^^jL = ad^^^^L + a(f){dH). 

The second Lie bracket on g = g © K, being an extension of (12. ip . is defined by 

[(a,a),(6,/5)]^ :=([a,6] R,^R{(i,b)) a,h E Q a,/9GK, (2.48) 

where 

ujR{a,h) =uj{Ra,h) +uj{a,Rh), (2.49) 

and Co" is a two-cocycle from (I2.4ip . Clearly, (I2.48P is a well-defined Lie bracket on g if and 
only if i? : g ^ is a classical _R-matrix and (I2.49p is a two-cocycle with respect to the second 
Lie bracket on g (12. ip induced by R. 

Proposition 2.13 A sufficient condition on R for (I2.49P to he a two-cocycle with respect to 
(12. ip is the Yang-Baxter equation (12. 3p . 

Proof. Skew-symmetry of (I2.49P is obvious. Hence, it suffices to verify the cyclic condition 
for the two-cocycles. Thus, from the Yang-Baxter equation (12. 3 p we have that 

cj/j ([a, 6]^ , c) + c.p. = u{R [a, fe]^ , c) + ([a, fe]^ , Rc) + c.p. 

= ijj ( \Ra, R]S\ , c) + auj ( [a, 6] , c) + ( [i?a, 6] , Rc) + ( [a, , Rc) + c.p. 
= uj{[Ra, Rb] , c) + uj{[Rb, c] , Ra) + uj {[c, Ra] , Rb) + c.p. = 0, 

where the last two equalities hold because a; is a two-cocycle. □ 
Hence, the second Lie- Poisson bracket on C°°(g) at a point (L, a) G g has the form 

{H, F}^ (L) := ((L, «), [(rfF, 0), (dH, 0)]^)^ 

= (L, [dF, dH]j^)^ + aujR {dF, dH) , (2.50) 

where H,F E C°°(g) and a G K. When a two-cocycle is given by (I2.44p . then the associated 
Poisson tensor ttr such that {H, F}^ = {dF, njidH)^ has the form 

TTndH = adJjrf^L + R*^!^^^ = ^'^*RdHL + i?*ad2^L + a(f){RdH) + aR*(f){dH). (2.51) 

Notice that the higher-order Poisson tensors from Sections 12.61 and | 2. 71 do not survive the pro- 
cedure of central extension. We have the following straightforward extension of Theorem 12.71 

Theorem 2.14 The Casimir functional Cn of (I2.47p . i.e., Cn G C°°(g) such that 

irdCn = ad^c^L + a <p{dCn) = 0, (2.52) 

are in involution with respect to (12.500 and hence generate the following hierarchy of mutually 
commuting Hamiltonian evolution equations on q: 

Ltr, = T^RdCn = £id*jiac„L + a (j){RdCn). (2.53) 
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Consider an important special case of fl2.44l) . when </> is a derivation with respect to 
additional continuous space coordinate. Consider once more a Lie algebra g with a non- 
degenerate symmetric ad-invariant scalar product (12.271) defined by means of a trace form 
tr on g. Assume now that g depends in a nontrivial fashion on an additional continuous 
parameter ?/ G S^, which naturally generates the corresponding current operator algebra 

g = C^{S\g) (2.54) 

of smooth maps from to g. On the current algebra we define the following modified trace 
form Tr : ^ K, such that Tr(a) := Jgi tr(a) dy, where a G g. Then the scalar product reads 

(a, 6)g := Tr(a6) = / tr(a6) dy a,b eg. (2.55) 

Thus, assuming that the derivative with respect to y is a derivation of the Lie bracket in g, 
i.e., (I2.46P with (p = dy holds, we can define the so-called Maurer-Cartan two-cocycle 

u (a, b) = (a, dyb)~ = / ti{adyb) dy a,b e g. (2.56) 

In this case the Casimir functions of the natural Lie-Poisson bracket on the centrally 
extended Lie algebra satisfy (12.521) in the form of the so-called Novikov-Lax equation 

[dCn,L] + ady{dCn) =0. (2.57) 

This follows from the invariance of the scalar product, as in this case ad* = ad. The differen- 
tials of Casimirs dCn generate the following Lax hierarchy (12.531) 

Lt„ = [RdCr,, L]+a dy{RdCn) = URdCn, (2.58) 

where the Poisson tensor (I2.5ip takes the form 

TTRdH = [RdH, L] + R* [dH, L] + a dy{RdH) + a R*dy{dH). (2.59) 



2.9 Dirac reduction and homotopy formula 

It often happens that we need to restrict the dynamics under study to a submanifold defined 
via some constraints. In such question arises of whether and how one can reduce 

the Poisson tensors. 

Assume that the (linear) phase space U = Ui (B U2 is spanned by ui G Ui and U2 G W2, 
i.e. u = (ui,U2)^. We will only consider the simplest case of the Dirac reduction given by 
the constraint U2 = c, where c e U2 is an arbitrary constant. In many cases considering 
such constraint is sufficient. Besides, more complicated constraints can always be reduced by 
change of dependent variables to several constraints of the above type. 

The Hamiltonian system with the Poisson tensor before reduction has the form 

uA ^ /^7rii(Ui,U2) 7ri2(Ui,U2 
U2/, l7r2i(Ui,U2) 7r22(Ui,U2 
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where we assume that 1122 is nondegenerate and hence invertible. Taking the constraint U2 = c 
into consideration we find that 

= (U2)t|^ = vr2i(ui, c)- h vr22(ui, cj — . 

OU-x OU2 

Thus, we can express |^ in the terms of |^ and put it into 

/ \ / \ 5H ^ 5H 

(ui)t = 7r2i(ui, c)- h 7r22(ui, c)- — =: vr (ui)- — . 

OUl OU2 OUl 

Hence the reduced tensor takes the form 

7r''<='^(ui) = 7rii(ui, c) - 7ri2(ui, c) ■ [7r22(ui, c)]"^ ■ 7r2i(ui, c). (2.60) 
Lemma 2.15 The operator f l2.60p is a Poisson operator on the affine space Ui © c. 

The skew-symmetry of f l2.60p is obvious. However the proof of the Jacobi identity for 
fl2.60l) consists of tedious yet rather straightforward calculations [22] and we will omit it. In 
the finite-dimensional case the said proof is much simpler (see for example [31]). 

Notice that if the inner product is not ad-invariant or we use the central extension pro- 
cedure, then in general we do not know the explicit form of the Casimir functions (like 
(12.291) for example) for the natural Lie Poisson bracket (or the extension thereof). In such 
a case one has to look for the annihilators dC of the Lie-Poisson tensor by directly solving 
the equation ndC = 0. With dC in hand, one can try to reconstruct the Casimir func- 
tions C. The Poincare lemma says that if the phase space L( is linear or of the star shape 
(Vn G U {eu :0^e^l}cW), then each closed /c-form is exact. In particular, when U satis- 
fies the condition of the Poincare lemma, we can reconstruct the Casimir functions C G C°^{q*) 
from their differentials dC G Q using the homotopy formula [61] 

Civ) 



[ {dC{er]),r])de r/ G 0*. (2.61) 
Jo 



Nevertheless, even when (12.611) is not applicable, the functions C can be often reconstructed 
through explicit integrations. 



2.10 Lax hierarchies from loop algebras 

Let g be a Lie algebra over K with the Lie bracket [-,■]. In contrast with Section [^721 we do 
not assume here any additional structures on g. 

Definition 2.16 The loop algebra over q is the algebra := g[A,A~^] of formal Laurent 
series in the parameter A G IK with the coefficients from g. 

It is easily seen that thanks to bi-linearity the operation in the former algebra extends to 
the loop algebra. Thus, in our case we can readily extend the Lie bracket [■, ■] to the loop 
algebra by setting 



[aA'", 6A"] = [a, 6] A'^^" a,beg m,neZ. 



(2.62) 
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There are two natural decompositions of into the sum of Lie subalgebras, i.e. 

0" = si©fl^ = I Vm,A^ I © I Vm,A^ I (2.63) 




for k = and 1. Thus for A; = and k = 1 we have well-defined classical /^-matrices (12.121) 

R = P+-l- (2-64) 

The transformation A t-^ A~^ maps the case oi k = into that of A; = 1, and vice versa. For 
this reason in what follows we restrict ourselves to considering the case of /c = only, and 
hence P+ and P_ will stand for projections onto nonnegative and negative powers of A. In 
fact, we have an infinite family of classical i?-matrices 

Rn = Ry' n e Z (2.65) 

and the corresponding new Lie brackets on g read 

[a, h]^^ := [Rna, b] + [a, , a, 6 G g. (2.66) 

The i?-matrices fl2.65p are well-defined since A*^ is an intertwining operator and Proposition l2.3l 
holds. 

Let L be an element of g. We have the following Lax hierarchy 

= [{X^L)^ ,L]=- [{X^L)^ ,L] neZ. (2.67) 

The commutativity of the flows (I2.67P for different n follows from the fact that (12.641) satisfy 
the Yang-Baxter equation (12. 3p and X„(L) = A"L are invariants (12.51) such that {X„i{L))^ = 
[RXm{L) , Xn{L)] still holds. The details of computations are parallel to those from the proof 
of Proposition 12. 4[ 

Without concentrating the specific properties of g, we can investigate the general form of 
appropriate Lax operators from g'^, i.e., the operators L E that generate self-consistent 
evolution equations on g'^ from the Lax equations (I2.67p . This means that the maximal 
and minimal orders in A of right- and left-hand sides of (12.670 have to coincide. Consider a 
bounded Lax operator L G g'^ of the form 

L = unX^ + un^iX^^^ + ... + Ui^m.X^""' + u^m^-"' iV, m G Z, (2.68) 

where ^ — m and Ui G g. Then a straightforward analysis shows that (I2.68P yields consis- 
tent equations (I2.67P if mat is a nonzero time-independent element of g, i.e., {uj^)t^ = 0. The 
specific properties of g might lead to further restrictions on (I2.68P . 

Theorem 12.71 can be applied for the construction of Hamiltonian hierarchies on the dual 
algebra to g^. We understand the dual algebra as g*^ = g*|A, A~^]. However, in contrast with 
the previously considered algebras, in the case of a loop algebra we have an infinite family of 
i?-matrices (12.651) with respective Lie-Poisson brackets (12.211) and related Poisson tensors of 
the form TTji^dH = ad^ ^^7] + R1ad*^jjri, where rj G Q*^ and dH G Q^. 

The coadjoint action ad* is defined with respect to the Lie bracket (12.621) on g'*'. All the 
above Lie-Poisson brackets are mutually compatible, which follows from the fact that the 
sum of intertwining operators also is an intertwining operator. Besides, if C„ G C°° {d*^) 
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is a Casimir function of the natural Lie-Poisson bracket fl2.19l) . i.e. a.d*^^^'!] = 0, then A^C„ 
is also a Casimir function. Hence, in the loop algebra case the Casimir functions generate 
multi-Hamiltonian hierarchies of mutually commuting vector fields on g*^ 

where dCn+i = X^dCn for / G Z. 

If we have a trace form on g, given by a linear map tr : g — *■ K such that the form in 
question is non-degenerate and symmetric, then this form can be extended to the loop algebra 
g^ by the formula 

IV(a) = tr(res a) a G g^, (2.69) 

where res a^A* = a_i. In fact, the choice of residue in (12.691) is a matter of convention, and 
one can choose the coefficient of an arbitrary order to get a proper definition of a trace. This 
is in contrast with the trace form (13.21) in the case of Poisson algebras. Nondegeneracy and 
symmetry of (12.691) are preserved. Moreover, if tr defines an ad-invariant metric on g, then 
this is also true for (I2.69p . and one can make an identification ad* = ad. 

Consider now the central extension approach for the loop algebras. Assume for simplicity 
that on g, and therefore on q^, we have a nondegenerate inner product, and hence g*^ = g'^. 
Let g^ be an extension of g^ with the Lie bracket (12.411) defined by a two-cocycle uj. The 
extended natural Lie-Poisson bracket has the form 

{H, F} (L) := (L, [dF, dH])^, + a uj {dF, dH) . (2.70) 

In contrast with the previous case, we have here an infinite family of /^-matrices (12.651) 
inducing an associated infinite family of new Lie-Poisson brackets on C°° (g^) , 

{H, F}„ (L) := (L, [dF, dH]J^ + a (rfF, dH) , (2.71) 

where 

ur^ (a, b) := UJ (i?„a, b) + uj (a, Rnb) a,b E g^. (2.72) 

All quantities (I2.72p are two-cocycles of the respective Lie brackets (I2.66p . For n = accord- 
ing to Proposition 12.131 this follows from the fact that R satisfies the classical Yang-Baxter 
equation, in the remaining cases one must additionally use the fact that A" are intertwining 
operators. It is important to stress that all Poisson brackets (I2.7ip are pairwise compatible, 
which follows from the fact that the linear sum of intertwining operators is an intertwining 
operator. If the two-cocycle u is given in the form (12.440 . then the Casimirs of the extended 
natural Lie-Poisson bracket (12.701) satisfy 

TTdCn = ad*ac„L + a <p{dCn) = 0, (2.73) 
and generate multi-Hamiltonian Lax hierarchy 

Lt„ = ad;^dC„-^ + " 4>{RdCn) = ... = 7lR_^dCn+l = TfRdCn = TlR^dCn-l = (2.74) 

where L E Q^, with the restriction that dCn+i = X^''dCn for I G Z. The respective Poisson 
tensors associated with the brackets (I2.7ip are given by 

TTR^dH = ad^^^^L + i^^ad^^L + a (f){RidH) + a Ri(j){dH). (2.75) 
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We readily see that R* = —A'/?. To find dCn, we can assume that 

dCn = X'^dCo = ao + aiX^ + 02-^"^ + . . . n e Z, 

and solve fl2.73|] recursively for the coefficients Oj. Then the functions C„ can be reconstructed 
using the homotopy formula fl2.6ip . 

Consider now a particular case of the Maurer-Cartan two-cocycle f l2.44p . with (p = dx and 
the loop algebra over (12.541) . Assuming that ad* = ad, the Lax hierarchy fl2.74p takes the 
form 

Lt^=[idCn)+,L]+adx{dCr,)+ = ...=Trn^dCn-i = ... • (2.76) 

Then, analyzing (12.760 we find that L G 0'*' of the form (12.680 yields self-consistent equations 
for ttT^OifA^^— 1 and Mat is nonzero and time- independent (except for = — 1) and 
m ^ 0. The Hamiltonian structures for (I2.76P are given by the Poisson tensors (I2.75P taking 
the form 



7iR,dH 



{X^dH) L -X^[dH,L]^ + adx{X^dH)+-aX^dx{dH)+. (2.77) 



The Poisson tensors (12.770 form a proper subspace of g'*' with respect to (12.680 . with the above 
restrictions, if ^ / ^ — m for A^ ^ and if ^ / ^ — m for A^ = — 1. Thus, there always exist 
at least two Poisson tensors ttr^ for which the procedure of Dirac reduction is not required. 
Recall that this analysis disregards specific properties of the Lie algebra q. Thus, if L is further 
constrained according to these specific properties the Dirac reduction might be yet required. 

Finally, let us note that in contrast with the previously considered algebras, where the 
central extension has lead to (2-1-1 )-dimensional systems, in the case of loop algebras the 
central extension is necessary for the construction of (l-l-l)-dimensional integrable continuous 
field systems. The reader will find the examples of this construction in the next sections. 

A natural choice are loop algebras defined over finite- dimensional semi-simple Lie algebras. 
In such a case the Killing form gives us symmetric, nondegenerate and also ad-invariant inner 
product. Thus, taking into consideration the central extension procedure with the Maurer- 
Cartan two-cocycle, the above choice leads to the constructions of a wide class of (1-1-1)- 
dimensional integrable continuous systems. The simplest case is presented in Section 14.31 



3 Integrable dispersionless systems 

The theory of integrable dispersionless or equivalently integrable hydro dynamic-type systems, 
i.e., the quasi-linear systems of first-order partial differential equations, belongs to the most 
recent ones and has been systematically developed from the 1980 's. Signiffcant progress was 
achieved after S. Tsarev discovered a technique called the generalized hodograph method 
that permits to find solutions using quadratures, see also [2H]- 

The study of the Poisson structures of dispersionless systems was initiated by B. Dubrovin 
and S. Novikov [25]. They established a remarkable result that the Poisson tensors of hydro- 
dynamic type can be generated by contravariant nondegenerate fiat Riemannian metrics. The 
natural geometric setting for the associated bi-Hamiltonian structures (Poisson pencils) is the 
theory of Frobenius manifolds based on the geometry of pencils of contravariant Riemannian 
metrics [21] • The Frobenius manifolds were introduced by B. Dubrovin as a coordinate-free 
form of the associativity equations, appearing in the context of deformations of 2-dimensional 
topological field theories (TFT), the so-called WDVV equations, that can be identified with 
(a class of) hydro dynamic-type systems. 
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3.1 Poisson algebras of Laurent series 

Consider the algebra of 'formal' Laurent series in p G C* about oo, 

Ki=k—r ) \i<k—T ) 

where Ui are smooth functions of continuous variable x G fi, i.e., A consists of polynomial 
functions in p and with finite highest orders, where f2 = if we assume these functions 
to be periodic in x or ^2 = M if these functions belong to the Schwartz space [ui and all their 
derivatives tend rapidly to zero when x approaches ±oo). We can introduce the Lie algebra 
structure on A in infinitely many ways using a family of Poisson brackets 

y .= pr r f g e A (3 1) 

' \dp dx dx dp J 

that generalize the well-known canonical Poisson bracket (the case r = 0). 

The trace form in the algebra A with fixed Poisson bracket (13.1 p for some r is defined in 
the following fashion: 

Tr/ := - / resoo (p~7) dx feA. (3.2) 

Here res is the standard residue at p = oo such that resoo-^^ = — M-i for L = '^^Uip^ 

Proposition 3.1 The scalar product (12.271) defined by means of the trace form (13. 2p is sym- 
metric, nondegenerate and a.d-invariant. 

Proof. The nondegeneracy and symmetry are obvious. Let 7 be a closed curve encircling 
once an infinity point on the extended complex plane. Then 

Tr{/,fi'}^ = - / Tes^{dpfd^g)dx + / ieSoo{d^fdpg) dx 
Jq Jn 

= / / ^^pdxg) dpdx + j j {d^fdpg)dpdx = 0, 

where the latter equality follows from integrations by parts with respect to p and x. Hence, 
the ad-invariance follows by Lemma (12.81) . □ 



Obviously, (13.10 is a derivation for the multiplication in A and hence we can further apply 
the scheme from Section 2.5. Fixing r we fix Poisson algebra and we are able to construct 
i?-matrices following from the decomposition of A. Simple inspection shows that A^k-r and 
A<:k-r are Lie subalgebras of A only in the following cases: 

• if r = for /c = 0; 

• if r e Z for A; = 1,2; 



• if r = 2 for /c = 



3. 
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Thus, fixing r we fix the Lie algebra structure with k numbering the i?-matrices given in the 
form R = P^k-r — \- Hence, we have multi-Hamiltonian Lax hierarchies ()2.33p of the form 



Lt„ = I (Iv^) , L^^ = TTodH^ = T^idHn-i n = 1, 2, . . . , (3.3) 

generated by fractional powers of infinite-field Lax functions L ^ A given in the form 

L = unp^ + + um-2p'''^ + ... iV 7^ 0, (3.4) 

where = l,UAr-i = for /c = and mat = 1 for /c = 1. A simple analysis of (13.31) shows 
that (13. 4p are appropriate Lax functions except for the case /c = 3, which is excluded from 
further considerations. This means that all the Lax functions under study have pole or zero 
at infinity. 

We find that R* = \- P^2r-k- The differentials of a given functional H e (A) of M 
have the form dH = Yli^^^-^"^ ^P^~^~^- Hence, the related Poisson tensors for (13.31) are given 
by (EBl 

{(LW)^,_,,L}^-Ln{rfi/,i^},)^2.-.' (3-5) 
where g = 0, 1, . . ., with the following Hamiltonians (I2.29P 

Hn = ^Tr (L^+i) n ^ -N. (3.6) 

We still have to check whether the above Lax functions span proper subspaces with re- 
spect to the above Poisson operators (13.51) . of the full Poisson algebras. We will restrict our 
considerations to linear in = 0) and quadratic {n = 1) Poisson tensors, as they are obvious 
enough to define bi-Hamiltonian structures. Besides, in all nontrivial cases, the Lax functions 
do not span proper subspaces w.r.t. the Poisson tensors for n ^ 2. 

For k = the above Lax functions always span the proper subspace w.r.t. the linear 
Poisson tensor, but for k = 1,2 this is the case only if ^ 2r — 2k + 1, otherwise the 
Dirac reduction is required. For the quadratic Poisson tensors the Dirac reduction is always 
necessary. The reduced quadratic Poisson tensor for /c = r = 0, 1, 2 is given by 

TTl^'^dH = {(LdH)^, ,L}^-L {{dH, L},)^^, + 1 {L, d'hes^ {dH, L},}^ , (3.7) 

see Section [23] and Lemma (I2.15p . The reduced Poisson tensors (13. 7p are always local as the 
residue from the last term is always a total derivative. 

The dispersionless systems with the Lax representations of the form (13. 3p where r = 
with k = 0,1 (the canonical Poisson bracket) related to the dispersionless KP hierarchy and 
the dispersionless modified KP hierarchy and the case /c = r = 1 of the dispersionless Toda 
hierarchy, together with their finite-field reductions, were considered in many papers (see for 
example [MlESllllllSlinnilHniElEZlIISland more recent pSl [IHl ED)- The theory with the 
Poisson bracket (13. ip for arbitrary integer r, from the point of view of classical /^-matrices, 
was considered for the first time in [8] and further developed in [121 EE] . 

The decomposition of the algebra A into Lie subalgebras is preserved under the trans- 
formation p I— > (the case k,r goes to k' = 3 — k,r' = 2 — r), but the Laurent series at 
oo (13. 4p transform into Laurent series at zero. This fact suggests more analytic approach to 
the construction of dispersionless systems, see [l6l [36] • For the theory of meromorphic Lax 
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representations (I3.3p of dispersionless systems see [SB]. In [IH] Krichever introduced the so- 
called universal Whitham hierarchies by means of the moduli spaces of Riemann surfaces of 
all genera. Other classes of reductions yielding nonstandard integrable dispersionless systems 
can be found for instance in [351 E21 EE] . For a more detailed discussion of the Lax represen- 
tations (13.31) for noncanonical Poisson brackets, their reductions and Hamiltonian structures 
as well as several examples of finite-dimensional reductions see [12] and [SB] . 

We know that the quantization of Poisson algebras for r = (and for the equivalent case 
r = 2) gives the algebra of pseudo-differential operators and leads to the construction of 
field soliton systems, while the quantization of the case for r = 1 gives the algebra of shift 
operators and leads to the lattice soliton systems. However, the class of reductions yielding 
construction of dispersionless systems is much wider then the class of corresponding reductions 
yielding systems with dispersion. Besides, the issue of quantization of Poisson algebras for 
r 7^ 0, 1,2 that would lead to the construction of dispersive integrable systems is still open, 
see [HH [HI [20], [HZ] and references therein. 

Below we present two examples of hydrodynamic chains (infinite-field systems) with their 
bi-Hamiltonian structures. For the classifications of hydrodynamical chains and related hy- 
drodynamic Poisson tensors see [Ml [SH [SSj • 

Example 3.2 Let us consider the infinite-field Lax function (13. 4p for N = 1: L = p + 
^^^qMjP"*"^. Then we find the first nontrivial hydrodynamic chain from hierarchy (13.31) . the 
well-known Benney moment chain, 

-f/ta = ^ I {L^)^q 5 = nodH2 = TiidHi {ui)t2 = (ui+i)^ + mi_i(uo)x, 

where {L'^)^o = + 2uq. The hi- Hamiltonian structure is given by the Poisson tensors 
T^o = jd^Ui+j^i + iui+j^idx 

-k'I = (z + \)Uij^jd^ + (j + l)<9a;^ij+j + + l)3Ui_xd^Uj_x 

i-i 

k=0 

which are obtained from (13. 5p for q = and (13.70 . The respective Hamiltonians are (13. 6p 

Hi = I- I ui dx H2 = ]r 1 [ul + U2) dx. 

Example 3.3 The case of k = r = 1 with the Lax function (13. 4p for N = 1 of the form 
L = p + Yli^o'^iP"^ ■ /irsi hydrodynamic chain from (13.30 has the following form: 

Lt^ = {(L)^o,-^}i = TTodHi = TiidHo <^=^ {ui)t^ = {ui+i)^ + iUi{uo)x- 

The bi-Hamiltonian structure is given by the Poisson tensors 

"^0 ~ j^xUi+j + iui^jdx 

i 

'^1 = [(j - k)ukdxUi+j_k + - k)ui+j^kdxUk] + i{j + l)uidxUj 

k=0 

-\- {i + l)dxUi+j+i + {i + l)ui+j+idx 
as well as the Hamiltonians Hq = f^u^dx and Hi = (ui + |mq) dx. 
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The central extension approach from Section 12.81 with the Mauren-Cartan two-cocycle 
yields the following Lax hierarchy fl2.58p of (2+l)-dimensional hydrodynamic systems 

Lu = {{dCn)^k-r^L~(^(l}r = ^odHn 71 = 1,2,..., (3.8) 

where the Poisson bracket takes the form 



{/jfi'lr- ^ ydp Qx dxdx)^\dqdy dy dq ) ^ ^ fi9^-^- 

The infinite- field Lax operators fl3.4p . with Ui depending additionally on ?/ e S^, are admissible 
with respect to fl3.8l) if ^ 1 — r. To construct evolution equations from fl3.8l) . we take 
dCr, = El =0 (^n-iP"' * and solve the Novikov-Lax equation (|2.57|), {dCn, L — aq}^ — for the 
auxiliary fields in terms of the fields from Lax function. The linear Poisson tensor (12.511) 
is given by 

TiQdH = {{dH)^^^^,L- aq]^- {{dH, L - a q}^)^^^_^ . 

The Lax functions (13.41) span the proper subspace w.r.t. the above linear Poisson tensor if 
^ 2r — 2fc + 1, otherwise the Dirac reduction is required. The construction of (2+1)- 
dimensional dispersionless systems by means of central extension procedure yielding the Lax 
hierarchy was presented in [13] , where one can also find a number of examples. 



Example 3.4 The case of k = r = 0. Consider L = +u. Then, for {dC3)^Q = + \up + 
jOid^^Uy we obtain the (2+l)-dimensional dKP equation 

3 3 

Lt^ = {((iC3)^o, L-aq}f^ = nodHs = ixidHi Ut,, = -uu^ + -a^d~^Uy, 

where the Poisson tensors are ttq = 28^^ and if a = ivi = d^u + ud^- The Hamiltonians are 
Hi = I I -V? dxdy = f f — (2u^ + 3a'^ud~'^Uyy) dxdy. 

Example 3.5 Consider the Lax operator L = p^^^ + up^~'^ + vp~'^ for k = 1 and r G Z, r 7^ 2. 
Then for {dC2-.r)^~r+i = p"^^^ + up^~^ we have 

Lt2-r = {{dC2-r)^l-r, L-aq}^ = 'KodH2-r = nidHi_r 

il\ _ (2 — r)vx + aUy 
V J ^ yrUxV + (1 — r)uv^ 

For r = 1 we get the (2+1)- dimensional dispersionless Toda equation with Hamiltonian struc- 
ture given by 

'''=\vdl oj and. If a = 0, vri=^ ^^^^ ^^^^^ 
with the Hamiltonians 

Hi = I I [ f + -Uj dxdy Hq = f f udxdy. 
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3.2 Universal hierarchy 

Let Q = Vect(S^) be the Lie algebra of (smooth) vector fields on the circle over the field 
K. The elements of Vect(§^) can be identified with smooth functions a{x) of spatial variable 
X G S^, with a Lie bracket in q of the form 

{a,b) := ab^ - ba^, (3.9) 

where a,b G g. Notice that 03.91] is a well-defined Lie bracket that does not satisfy the 
Lebniz rule 02.41) . Thus, we follow the scheme from Section P2. 101 and consider the loop algebra 
over Vect(S^), i.e., = Vect(§^)[A, A~^]. The commutator 03. 9p readily extends to g^. We 
already know that we have the decomposition of g^ into Lie subalgebras 02.631) . with respective 
classical i?-matrices generating the following Lax hierarchies 02.671) : 

L,„ = ((A"L)^,,L> A; = 0,1, (3.10) 

where L G g^ and n G Z. 

The hierarchy 03. 101) is the so-called universal hierarchy of hydrodynamic type, which is 
a subject of intensive research in recent years [55l [56]. The hierarchy (13.101) can be obtained 
as a quasi-classical limit of the coupled KdV equations of Antonowicz and Fordy [1]. In 
this fashion we can obtain multi-Hamiltonian structure for the universal hierarchy [30]. In 
fact, the multi-Hamiltonian structure of the coupled KdV equations using classical i?-matrix 
approach can be algebraically interpreted as a set of compatible Lie-Poisson structures on the 
dual space to the loop Virasoro algebra [32] . The Virasoro algebra is a central extension of the 
Lie algebra of vector fields on a circle Vect(S^) associated to the Gelfand-Fuchs two-cocycle 
(i.e., ([231 with = d^). 

Example 3.6 Consider (I3.10p for the infinite-field Lax functions given in the following (ap- 
propriate) form L = uo + miA~^ + M2A~^ + . . ., where uq = I for k = 0. One can observe that 
Lt„ = ((A"L)^^ , (A"L)^^) A~". Hence, the evolution equations from (13.101) take the form of 
(l+l)-dimensional hydrodynamic chains 

i-l+k 

{Ui)t„ = (1 - k){Ui)^ + ^ {Ui_j, Un+j) k = 0,l. 

j=i-k 

Assume that the dynamical fields in g'^ depends on an additional spatial variable y. Let 
us now consider a (2+l)-dimensional counterpart of (I3.10p . It reads 

Lt^ = {(An)^, ,L)+dy (A„)^, /c = 0, 1, (3.11) 

where An = a^A"" -|- Qn-iX"""^ + ■ ■ ■ satisfy 

{An,L) + {An)y = 0. (3.12) 

Notice that An = X^Aq. For a given L G g'^ one finds coefficients aj from An by solving 
03.12|) recursively. Notice that one cannot obtain 03. lip as a central extension of the universal 
hierarchy, as 02.56P is not a two-cocycle associated with Vect(S^). Commutativity of the 
equations from the hierarchy (13. lip can be proved by straightforward computation. On the 
other hand, integrability of the equations from (13. lip follows from the fact that (13. lip is a 
Lax sub-hierarchy of the centrally extended cotangent universal hierarchy |76j . 
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Example 3.7 Let L have the form L = X + u in the case of k = 0. Then solving fl3.12l) one 



finds that Aq = 1 + uX ^ + ^UyX ^ 



Hence, the first nontrivial (2+l)-dimensional 



hydrodynamic equation from the hierarchy (13.111) is f5^ 



't2 



\{A^)^,,L) + dy{A^) 



d,,. ^Uyy — UUy + UxO^ 



This hydrodynamic equation is equivalent to system of the form 

UU -Vy+ UVx - U^V =0 Vr, -Uy = 0, 



that has recently attracted considerable attention, see [53, [7^ . 

Example 3.8 The case of k = 1 for L = uX~^ . We have Aq = 1 — dy^u^X"^ + . . .. Hence 



't2 



{{A,)^^,L)+dy{A^) 



ut., = udy ^Uxx - Uxdy ^Ux. 



4 Integrable dispersive systems 

In the following section we apply the general formalism to several infinite dimensional Lie 
algebras in order to construct a vast family of integrable dispersive (continuous and discrete 
soliton) systems. 

Recently, the so-called integrable g-analogues of KP and Toda like hierarchies become of 
increasing interest, see [39l [331 HQl [3], EH] . Our approach presented in two following subsections 
includes the g-systems as a special case. Actually, we consider generalized algebras of shift 
and pseudo-differential operators that allows to construct in one scheme not only ordinary 
lattice and field systems, but in particular also their g-deformations. 

We also present (2 + l)-dimensional extensions of the Lax hierarchies following from the 
algebras of shift and pseudo-differential operators. In these cases the the quadratic Poisson 
tensor is not preserved. Nevertheless, the second Poisson tensor can be given by means of the 
so-called operand formalism [THl [HI [51] . Its construction, within classical i?-matrix approach, 
can be found in [T5] . 

In last subsection we present the application of classical i?-matrix formalism to the loop 
algebra s/(2, C)[A, A~^], which is the simplest case of infinite-dimensional Lie algebras of 
Kac-Moody. In fact the scheme of construction of infinite-dimensional systems from affine 
Kac-Moody Lie algebras is one of the most general and particular cases are closely related to 
the algebras of shift and pseudo-differential operators yielding soliton systems. For details we 
send the reader to the survey [23] • 

For the application of i?-matrix formalism to other algebras see for example [TQl [691 EHl [ZH] • 
A specially interesting class of algebras is related with the so called super-symmetric (SUSY) 
systems. The details the reader finds in [^Hl [HZl [IZl ISHl [SH] and in literature quoted there. 

4.1 Algebra of shift operators 

Consider the algebra of 'formal' shift operators 
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where Ui G and JF is an algebra of dynamical fields with values in K. The associative 
multiplication rule in g fl4.ip is defined by 

£™M = m G Z. (4.2) 

Proposition 4.1 The multiplication in (14. ip of the form fl4.2p zs associative if and only if 
E : T ^ T is an invertible endomorphism, i.e. 

E{uv) = E{u)E{v). 

The proof is straightforward. The Lie bracket in g is given by the commutator [A, B] = 
AB - BA, where A,B e Q. 

Let tr : g — s> K be a trace form, being a linear map, such that 

tr(A) := (free(A)) , (4.3) 

where free(A) := oq for A = Oji^* and (■) denotes a functional whose form depends on the 
realization of the algebra JF and the endomorphism E. We assume that (■) is such that 

{Ef) = (/) feJ' (4.4) 

holds. Notice that from this assumption it follows that (uEv) = {vE^^u), thus the adjoint of 
operator E is E'^ = E~^ . 

Proposition 4.2 The bilinear map defined as 

{A,B)^:=ir{AB) (4.5) 
is an inner product on g which is nondegenerate, symmetric and ad-invariant. 

Proof. The nondegeneracy of ( 14.50 is obvious. The symmetricity follows from the defini- 
tion by using (14. 4p . The ad-invariance is a consequence of the associativity of multiplication 
operation in g, see Lemma [2.81 □ 

Let us mention now a few standard realizations of the algebra of shift operators (14. ip : 

• The first one is given by the shift operators on a discrete lattice. In this case the 
dynamical functions are : Z — > K and E^u{n) := u{n + m). The form of the 
functional from (14.30 is (/) := Xlnez /("-)• 

• The second realization is given by the shift operators on a continuous lattice. Now 
Mj : f2 — > K are smooth functions of x, where f2 = §^ or f2 = M, if Ui are from the 
Schwartz space. In this case the shift operator is E"^u{x) := u{x + mh), where ^ is a 
parameter. In this case the functional is just integration, i.e., (/) := f{x)dx. 

• The third realization is for q-discrete functions : Kg ^ K, where Kg := U {0} for 
q^O. Here E'^m(x) := uiq'^x) and (/) := T.n&fil )■ 
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All the above realizations lead to the construction of lattice soliton systems: discrete, 
continuous and g-discrete, respectively. More realizations can be made by means of the 
discrete one-parameter groups of diffeomorphisms [9] or by means of jump operators on time 
scales [371 E]- Notice that in the continuous limit the algebra of shift operators on lattices 
gives the Poisson algebra for (13.11) with r = 1. Thus, the quasi-classical limit of lattice soliton 
systems are respective dispersionless systems. The situation in the g-discrete case is similar 

i- 

The subspaces Q^k-i and Q<k-i of (14. ip are Lie subalgebras only for k = 1 and k = 2 and 
the classical i?-matrices following from the decomposition of g are (12.121) R = P^k-i — \- Their 
adjoints with respect to the above inner product are given by R* = P^2-k — \i respectively. 
As a result, we have two Lax hierarchies: 

(L^) , l] = TlodHn = 7TldHn-l k = l,2, (4.6) 

of infinitely many mutually commuting systems. Let (14. 6 p be generated by powers of appro- 
priate Lax operators L E g oi the form 

L = un£^ + UN-i£^-^ + + UN-sS""-^ + (4.7) 

where ujv = 1 for A; = 1. 

The bi-Hamiltonian structure of the Lax hierarchies (14.60 is defined by the compatible (for 
fixed k) Poisson tensors given by the formulas (I2.34p . 

TTodH = [L, {dH)^k-i] + {[dH, L])<2_fc 

and 

TT.dH = 1 ( [L, {LdH + dHL)^,_,] + L {[dH, L])^,_, + {[dH, L])^,_, L) 

+ (2 - A;) [{E + 1){E- l)-i free {[dH, L\) , L] , (4.8) 

where the operation [E — is the formal inverse of {E — 1). The second Poisson tensor is 
a Dirac reduction of (12.380 as in this case R = R + {k — f-Po) satisfies the related condition. 

The differentials dH{L) of functional H{L) G C°^{q) for (gTD have the form dH = 
Yl'^2-k Sun respective Hamiltonians (I2.29P are 

HniL) = -^tr N ^ -n. 

The theory of lattice soliton systems of Toda type with Lax representations given by 
means of the shift operators was introduced for the first time by B. Kupershmidt in |47j . 
This class of systems, follows from the classical i?-matrix formalism applied to the algebra of 
shift operators, was investigated from this point of view in [10] and [57] . 

Example 4.3 Consider the case of k = 1 with (14.70 normalized as L = S + Yl^o UiS~\ The 
first chain from the Lax hierarchy (14. 6 p has the form 

Lt, = [{L)^^ , L] = TiodH^ = TiidHo ^ {ui)t, = {E - l)ui+i + ^^(1 - ^"0%, (4.9) 
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and its explicit bi-Hamiltonian structure is given by 
ttq^ = E^Ui+j - Ui+jE~' 

i 

= Yl [ukE^~'"ui+j-k - Ui+j^kE^-'uk + Ui {E^-^ - E~^) Uj\ 

+ Ui (l - E^~') Uj + E-^'+^n.+j+i - Ui+j+i^-^-^ 

together with the Hamiltonians Hq = {uq) and Hi = (ui + \u^ . 

Assume now that the dynamical fields depend on an additional variable ?/ G S^. Then 
after the central extension procedure with the Maurer-Cartan two-cocycle f l2.56p the (2+1)- 
dimensional Lax hierarchy takes the form f l2.58p . [T5] . 



Lt„ = [{dCn)^k-i ,L-ady] = nodCn k = 1,2, 

for the Lax operator (14 .yp with > 1, where dCn = X^i^o solutions of (I2.57p . 
[dCn, L — ady] =0. The Poisson tensor (12.590 is given by 

TTodH = [L- ady, {dH)^k~i\ + {[dH,L -ady])^^_^,, 

and the Dirac reduction is not required. 

Example 4.4 The case of k = 1. An example of finite field reduction. The Lax operator is 
given by L = £ + u + v£~^ . Then for {dCij^Q = S + u we have 

Lt, = [idCi);^o, L - ady] = ivodHi = TndHo ^ (^^^ " (^^(1 -E-^)u ") ' 

The respective Poisson tensors are 

f ady (E-l)v\ , ^ / Ev-vE~^ u(E-l)v 

^0= n TP-i\ n and, if a = 0, tti = , ' 

\^f(l — -6/ ) J \v [1 — E ^)u V [E — E 

The Hamiltonians are Hi = (^v + ^u"^) and Hq = (u). 

4.2 Algebra of (5-pseudo-differential operators 

Consider a generalized derivative in the algebra JF of dynamical fields with values in K given 
by a linear map IS.: T ^ T that satisfies the generalized Leibniz rule 

lS{uv) = A{u)v + E{u)A (4.10) 

for an algebra endomorphism E : IF ^ IF. If E' = 1, then A is an ordinary derivative. 
According to (I4.10p we define a generalized differential operator 

5u = A{u) + E{u)5. (4.11) 

In (14. lip and in all subsequent expressions A and E act only on the nearest A-smooth 
function on their right-hand side. As above, the expression (14. lip is a counterpart of an 
ordinary differential operator d such that du = u^ + ud. Using (14. lip we have 

5~^u = E~^u5~^ + 5~^A'^u5~^ = E^\5^^ + E-^A'^ud^^ + E~^A'^^u5~'^ + ..., 
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where := -AE'^. 

Hence, we can further define an algebra of 5-pseudo-differential operators 



[i^k J I i<k J 

where Ui G JF. The above algebra is noncommutative and associative. In fact we have the 
following result. 

Proposition 4.5 The algebra (I4.12p generated by the rule of the form (14.111) is associative 
if and only if A : T ^ T satisfies (14.101) and E : T ^ T is an algebra endomorphism. 

We omit the proof, which can be done by induction, and it suffices to consider the asso- 
ciativity condition, {AB)C = A{BC), only for A = 6\ B = bS^ and C = c, where a,b,c E T . 

In fact we will consider only two most important special cases of generalized derivatives 
and (5-pseudo-differential operators: 

• The first case is the ordinary derivative, i.e., (I4.10p with E = 1. Thus let !F consist of 
smooth functions : f2 — *■ IK of x, where = §^ or f2 = M (if Ui belong to the Schwartz 
space). Thus A = d^, and we let S := d. In this case q is the algebra of standard 
pseudo-differential operators with the trace form given by 

iiA = (lesA) = / lesAdx A = V^a^S* G q, resA := a_i. (4.13) 
Jn i 

• The second case is the generalized derivative given by a difference operator. Thus, let 

A = -(E-1) =^ 5 = -(S-l), 
/i /i 

where fi is constant, E is an algebra endomorphism and £" is a shift operator such 
that (14.21) holds, see the realizations from the previous section. In particular, in the 
continuous lattice case: 

A/(3;) = Z(^±^^l/M, Efix) = fix + h), fi = n; 

and in the g-deformed case: 

A/(x) = ^^MziM, Ef{x) = f{qx), /i = g-l. 
q-1 

Here the trace form is exactly the same as in the case of shift operators algebra from 
previous section (14.31) . Thus, we have to consider the restriction (14.41) again. Notice that 
(5-pseudo-differential operators can be represented uniquely by shift operators, with the 
convention that the (5-operators of negative orders are expanded into shift operators of 
negative orders as well. Note that in this case = —AE~^. 
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In the quasi-classical limit the algebra of standard pseudo-differential operators gives the 
Poisson algebra with canonical Poisson bracket fl3.1l) for r = 0. Thus, in dispersionless limit 
continuous soliton systems yields respective dispersionless systems. In the case when difference 
operator A is given on lattices, the continuous limit of the algebra of pseudo-5-differential 
operators is actually given by the algebra of standard pseudo-differential operators. Thus, 
in contrary to the previous case, the quasi-classical limit of discrete soliton systems yields 
continuous soliton systems. This situation is similar in the g-deformed cases |11] . 

The reason why we consider both cases, i.e., continuous and discrete, is the fact that they 
can be unified into a single consistent scheme using the so-called time scales [371 El 188] . 
This scheme also includes soliton systems with spatial variable belonging to the space being 
a composition of continuous and discrete intervals. In [HH] the trace formula on the algebra 
of pseudo-5-differential operators is unified to the form that covers all the above realizations 
as well the cases with nonconstant fi. 

Proposition 4.6 The inner product f l2.27p on g defined by means of traces f l4.13p and (14.31) 
in both cases is nondegenerate, symmetric and a.d-invariant. 

Proof. In the first case of ordinary derivative nondegeneracy is obvious, the symmetricity 
follows from integration parts and the assumption that integrals of total derivatives vanish. 
The ad-invariance is a consequence of Lemma 12.81 In the second case of difference operators 
the proof follows from Proposition 14.21 and the fact that 5-pseudo-differential operators can 
be expanded by means of shift operators. □ 

The Lie structure on (14.121) is given by the commutator [A, B] = AB—BA, where A,B ^ q. 
In general the subspaces g^k and g<fc are Lie subalgebras of g only for k = and k = 1. 
However, if E = 1 they are Lie subalgebras also for k = 2. The classical i?-matrices following 
from the decomposition of g are R = P^k — \- Hence, we have the following Lax hierarchies 
of commuting evolution equations 

= TldHn = TTldHn^l k = 0,l (4.14) 

and for k = 2 ii E = 1. In (14.141) L G 0. The infinite- field Lax operators generating (14.141) 
are given in the form 

L = + Uiv-i^"^"' + + njv-a^'^"' + . . . iV ^ 1, (4.15) 

where in general for k = the field is time-independent; if E = 1 then for k = the field 
mat-i is also time-independent and for k = 1 the same for the field u^. 

The explicit form of the differentials dH = 5~*~^7i with respect to general Lax operator 
(14.151) has to be such that (12.161) is valid. See the Example 14.81 In the case of = 1 we have 
7j = 1^. The Hamiltonians (12.291) are given by 

HniL) = -^tr N ^ -n. 

We will consider the Hamiltonian structures of (14.141) only for k = 0. Thus, for A; = we 
have R* = —R. Hence the linear Poisson tensor is given by (12.341) . 



(n \ 



7iodH=[{dH)^o,L]-{[dH,L]) 



(4.16) 
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The case of quadratic Poisson tensor is more complex. For the pseudo differential operators, 
when E = 1, the quadratic Poisson tensor is given by 



1 



nl^'^dH = {LdH)^Q L-L {dHL)^^ + — [9"^ (res [dH, L]),L] 



(4.17) 



after Dirac reduction applied to (12.381) . In the case of difference operators (-E 7^ 1), as the 
decomposition for = of the algebra of (5-pseudo-differential operators (I4.12p coincides with 
the decomposition of the algebra of purely shift operators (14.11) for /c = 1, the quadratic 
Poisson tensor is given by (14. 8p . see [23, EH] for further details. 

The pseudo-differential case {E = 1) of (I4.14p with k = and finite-field reductions, 
given by constraint Ui = for z < 0, is the well-known Gelfand-Dickey hierarchy [3l]. More 
details of the i?-matrix formalism applied to the algebra of pseudo-differential operators for 
the remaining values of k can be found in [13] or [7] . 

Example 4.7 Consider the infinite-field case k = of (I4.15P in the ordinary derivative case 
E = 1, i.e., L = d + ■ This is the case of the well-known KP hierarchy JH]/. 

Then we find the first nontrivial dispersive chain from (13. 3p . 



^t2 



where (L^ 



nodH2 = TTidHi 



{Ui)t2 = {Ui)2x + 2(Mj+l)x - 2 ^(-l)'' ( ^ j Ui-kiuo)kx 



k=l 



+ 2uo. The bi-Hamiltonian structure is given by linear 



E 

k=l 



i+j-k 



El 

fc=i 



-1^ 



Ui+j-kd^ 



and quadratic Poisson tensor 
i+i 



k=l 



k=l 



i + 1 

k 



Ui+j-k+id^ 



E E 

1=0 I k=i 

EEE(-i 

1=0 k=0 s=l 



J-l-1 
k 



uid^Ui+j-k-i-1 



k=l 



J-l-1 

s 



U^^k+l^^x^' 



that are obtained from (I4.16P and (I4.17p . respectively. The respective Hamiltonians are 



Hi 



Ui dx 



{ul + U2) dx. 



Example 4.8 The case ofk = 0. In this example we present the associated integrable systems 
and their Hamiltonian structures in the form that is valid in both continuous and discrete 
cases. In the first case, E = 1, fi = and A = — A"!" = d^. 
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Let the Lax operator be given in the form 

L = 6 + fiijjip + ijS~'^(p. (4.18) 
Then the first and the second flows from the Lax hierarchy (14 .6^ are 

and 



(4.19) 



(4.20) 



For the Lax operator (I4.18P the differential of a functional H such that (I2.16P is valid is 
given by 

,H = - iAt a) A" A - ^ ',„ A-M. 

(f dip %p \f J ipip + fiipA'if 

where A = — ^j^, o,nd A~^ is a formal inverse of A. The linear and quadratic Poisson 
tensors take the form Jg^ . 



_ l\ _ /-/iV'^ - 2?/^A-V A + 2^?/^v9 + 2V'A- V' 

"""-^ -1 oj' ~ 1^ -At + 2V9A-V _ 2<^A-V 

while the Hamiltonians are 

H2 = ^^/iVv'^ + ^/'V^ + V^^V + Ai^V'^A^ + /i^V^V 

/n particular, when E = 1, fi = and A = dx the above bi-Hamiltonian hierarchy is 
precisely the bi-Hamiltonian field soliton AKNS hierarchy J60^ . In this case the first nontrivial 
flow is the second one fl4.20p . i.e., the AKNS system. When A is the difference operator, we 
obtain in particular the lattice and the q-discrete counterparts of the AKNS hierarchy, 
where the first nontrivial flow is (I4.19P . 

In fact this example is more general and also includes more complex situations when /i is 
non-constant (time-independent) function on M, for the details we send the reader to f88^ . 

Let g be the algebra of pseudo-differential operators with the coefficients depending on 
independent spatial variables x and y. The central extension procedure with fl2.56p yields the 
following Lax hierarchy (l2.58p[T5] 



Lt„ = [idCn)^k , L - ady] = TiodCn, 

where k = 1,2 oi k = 3 and L = unQ^ + UN-id'^~^ -\-UN-2d^~^ + . . ., with mat = 1,un-i = 
for k = and only m^v = 1 if /c = 1. The Lax hierarchies are generated by dCn = ^n-id"'~^ 
solving (12.570 ■ [dCn, L — ady] = 0. The associated Poisson tensor (I2.59P is given by 

TiodH = [{dH)^o, L - ady] - {[dH, L - ady])^^^ , 

and requires no Dirac reduction. 
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Example 4.9 The case of k = with the Lax operator of the form L = d"^ + u. Then, for 
{dC3)^o = + + I + ad'^Uy) we obtain the (2+1)- dimensional KP equation 

13 3 

= [(rfC3)jjo, L - ady] = ixodH-i = TXidHi <^=^ = -u^xx + -uux + -a'^d'^Uy, 

where the Poisson tensors are ttq = 2dx and, if a = 0, tti = ^d^+dxU+udx- The Hamiltonians 
are 

Hi = -v^ dxdy = — (2u^ + uu^^ + 3a'^ud~'^Uyy^ dxdy. 

More examples of (2 + l)-dimensional field systems can be found in fTEl [Mj . where one 
can find also systems that are purely (2 + l)-dimensional phenomena. 



4.3 s/(2,C) loop algebra 

We will follow the scheme from Section I2.10[ Consider a loop algebra over the classical Lie 
algebra g = s/(2,C) of traceless 2x2 nonsingular matrices, i.e., = s/(2, C)|A, A~^], with 
coefficients being smooth dynamical functions of variable x E Q. The space Ui : Q ^ K are 
smooth functions of x, where f2 = §^ if we assume these functions to be periodic in x or ^2 = M 
if these functions belong to the Schwartz space. The commutator defines the Lie bracket in 
s/(2,C) and readily extends to since (12.621) . 

We already know that there are two natural decompositions of into Lie subalgebras 
(12.631) . and we consider only the one for k = yielding the classical i?- matrix i? = P+ — |, 
with P+ being the projection onto the nonnegative powers of A. 

The trace form on g^ is (12.691) . Tr(a) = restr(a) dx, where res Yli '^j-^* = '^-i tr is the 
standard trace of matrices. As the matrices in s/(2, C) are nonsingular, the trace tr defines 
nondegenerate inner product which is also symmetric and ad- invariant. These properties 
extend to g^, and hence g*^ = g^ and ad* = ad. 

Applying the central extension procedure with the Maurer-Cartan two-cocycle (I2.56P (x = 
y) yields the Lax hierarchy given by (I2.74p (we take a = 1), i.e., 

Lt^ = [{dCn)+ ,L-dx]=...= TTidCn-l = . . . (4.21) 

for the Lax operators 

L = unX^ + W7v-iA^~^ + . . . + Wi^^A^-™ + w_„A-™ N ^ -1, (4.22) 

where Ui G s/(2,C) and is a constant matrix. The Lax hierarchy (I4.2ip is generated by 
dCn = X^dCo such that dCo = Ylilo where G s/(2, C), satisfies [dCo, L — dx] = 0. 
For a given functional H G C°°g^ of (I4.22p its differential has the form 

- -A"*"^ + A"'"^ + . . . + A"^, 



where 



6H 



1 SH SH 

, 'm' 1 for = ( 
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The multi-Hamiltonian structure for fl4.21l) is given by fl2.77l) 



TiidH 



{X'dH) L-d^ -X'[dH,L-d.,]^. (4.23) 



For the general Lax operators (14.221) . ii N ^ I ^ — m for N ^ Q and ^ / ^ — m for = — 1 
then the Dirac reduction of (14.231) is not required. 

Example 4.10 Consider the Lax operator of the form 



i 0\ , (Q q 



Then we find that 



Hence, from (14.211) we obtain the Ahlowitz-Kaup-Newell-Segur (AKNS) hierarchy JI]/ 

Vti VW ' V)t2 \-^r^^ + ir^qJ ' \rj^^ V-i^^^^ + l'"^'"^ 
The bi-Hamiltonian structure is given by the Poisson tensors (14.231) 

2qd^^q — '2qd~^r 





-2A 




\2i 


j 





^° \2t J \d^-2rd^^q 2rd~^r 

with the hierarchy of Hamiltonians 

Hi= [ ^qxr dx, H2 = [ ^ [q^r"^ - g^^r) dx, H3 = [ ^ {^Qxt"^ - qxxxr) dx, ... . 

Notice that the need for the Dirac reduction for the Poisson tensor ttJ*^'^ follows from the fact 
that (I4.24P is not the most general Lax operator (I4.22p from with N = 1 and m = 0. 

The reduction q = ip ^ r = ip* of the above AKNS hierarchy yields the non-linear Schrddinger 
hierarchy, while the reductions q = iu,r = i and q = r = iv give rise to the KdV and the 
mKdV hierarchies respectively. 
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